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ABSTRACT 

The ideal helical kink instability of a force-free coronal magnetic flux rope, anchored in the photosphere, is 
studied as a model for solar eruptions. Using the flux rope model of Tito v & D emoulin ( 1999) as the initial 
condition in MHD simulations, both the development of helical shape and the rise profile of a confined (or 
failed) filament eruption (on 2002 May 27) are reproduced in very good agreement with the observations. By 
modifying the model such that the magnetic field decreases more rapidly with height above the flux rope, a 
full (or ejective) eruption of the rope is obtained in very good agreement with the developing helical shape and 
the exponential-to-linear rise profile of a fast coronal mass ejection (CME) (on 2001 May 15). This confirms 
that the helical kink instability of a twisted magnetic flux rope can be the mechanism of the initiation and the 
initial driver of solar eruptions. The agreement of the simulations with properties that are characteristic of many 
eruptions suggests that they are often triggered by the kink instability. The decrease of the overlying field with 
height is a main factor in deciding whether the instability leads to a confined event or to a CME. 

Subject headings: Instabilities - MHD - Sun: corona - Sun: flares - Sun: coronal mass ejections (CMEs) 



1. EVITRODUCTION 

' Large-scale solar eruptions occur as flares, filament (or 
prominence) eruptions, and coronal mass ejections (CMEs). 
Despite their different observational appearance, it is believed 
that these phenomena are manifestations of the same physical 
processes, which involve the disruption of the coronal mag- 
netic field. Indeed, in the largest eruptions (eruptive flares) 
usually all three phenomena are observed. The theory of the 
main phase of such events, referr ed to a s the "standard model" 
of eruptive flares (e.g., Shibatal fl 999ft . is quite well estab- 
lished. However, their initiation as well as the mechanism of 
upward acceleration are still unclear A variety of theoretical 
models have been proposed to explain the impulsive onset and 
initial evolution of solar eruptions (see, e.g., Forbes 2000). 

Here we focus on a flux rope instability model. This is mo- 
tivated by the observation that solar eruptions often show the 
phenomenology of a loop-shaped magnetic flux system with 
fixed footpoints at the coronal base and signatures of magnetic 
twist. Furthermore, erupting filaments very often develop a 
clearly helical axis shape in the course of the eruption, which 
is the characteristic property of the helical kink instability of a 
twisted magnetic flux rope. The instability occurs if the twist, 
a measure of the winding of the field lines about the flux rope 
axis, exceeds a critical value (Hood & Priest 1981). 

In coronal applications, the simplifying assumption of 
straight, cylindrically symmetric flux ropes has nearly always 
been used so far. Only very recently, iTorok. Kliem. & Titovl 
( 12004 hereafter Paper I) performed the first detailed study 
of the kink instability of an arched flux rope, line-tied to the 
photosphere, using the analytical model of a force-free coro- 
nal flux rope developed by Titov & Demoulin ( 1999, hereafter 
TD) as the initial condition in 3D ideal MHD simulations. 
They have shown that this model relaxes to a numerical equi- 
librium very close to the analytical expressions in the case of 
subcritical twist and that the helical kink instability develops 
for supercritical twist (see also Fan & Gibson 2003, 2004). 
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In the course of the instability, a helical current sheet, 
wrapped around the kinking and rising flux rope where it 
pushes into the surrounding field, and a vertical current sheet 
below the rope (which has no counterpart in the cylindrically 
symmetric case) are formed. A vertical current sheet below 
rising unstable magnetic flux is the central element in the stan- 
dard model of eruptive solar flares. Further essential features 
of solar eruptions could be reproduced in the simulations, as 
for example the formation o f transient soft X-ray sigmoids 
dKliem. Titov. & Torokl2004l Paper II). However, a full erup- 
tion of the configuration has not yet been obtained; the flux 
rope reached an elevation of only about twice its initial height. 

Here we present further developments of these simulations 
to substantiate our suggestion in Papers I and II that the kink 
instability of a coronal magnetic flux rope is a possible trig- 
ger mechanism of solar eruptions. The instability was first 
suggested as the trigger of (confined and ejective) prominence 
eruptions by Sakurai ( 1976), but has recently been generally 
regar ded as a po ssible explanation only for confined events 
(e.g..l Gerrard & Hood.2003.) . The new simulations show that 
the instability can also trigger full eruptions. 

2. NUMERICAL MODEL 

We integrate the compressible ideal MHD equations using 
the simplifying assumptions of vanishing plasma-beta, (3 = 0, 
and vanishing gravity, which are identical to Eqs. (2-5) in Pa- 
per I. Setting /3 = is usually a very good approximation in the 
lower and middle corona of active regions, the source region 
of most eruptions, where estimates yield (3^ 10"^... 10"^. 
Both the pressure gradient force and the gravity force influ- 
ence the rise characteristics of the unstable flux rope and the 
energy partition in the development of the instability. How- 
ever, the basic characteristics of the instability are well de- 
scribed by the equations used whenever the Lorentz force 
dominates, as is the case in the initial phase of solar erup- 
tions. Magnetic reconnection occurs in the simulations due to 
numerical diffusion if current sheets steepen sufficiently. 

As in Papers I and II, we use the approximate analytical 
force-free equilibrium of an arched, line-tied, and twisted flux 
rope by TD as the initial condition for the magnetic field. 
The flux rope is modelled by the upper section of a toroidal 
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Fig. 1. — Left: TRACE 195 A images of the confined filament eruption 
on 2002 May 27. Right: magnetic field lines outlining the core of the kink- 
unstable flux rope (with start points in the bottom plane at circles of radius 
b/3) at ? = 0, 24, and 37. The central part of the box (a volume of size 4') is 
shown, and the magnetogram, B;(ji:,j',0,f), is included. 



ring current, partly submerged below the photosphere, whose 
Lorentz self-force is balanced by a pair of fictitious subpho- 
tospheric magnetic charges. A fictitious subphotospheric line 
current at the toroidal symmetry axis is included to achieve a 
finite twist everywhere in the system. See TD for a detailed 
description of the model. 

The initial density distribution can be freely specified; we 
choose it such that the Alfven velocity in the volume sur- 
rounding the flux rope decreases slowly with height: po oc 
|Bo(x)p''^ (see Fig.|3]below). The system is at rest at f = 0, 
except for a small upward velocity perturbation, which is lo- 
calized at the flux rope apex in a sphere of radius equal to the 
minor radius, b, of the rope. Lengths, velocities, and times are 
normalized, respectively, by the initial flux rope apex height, 
ho, the initial Alfven velocity at the apex, vao, and the corre- 
sponding Alfven time, ta = /iq/vao- 

3. SIMULATION OF A CONFINED ERUPTION 

The eruption of an active region filament (on 2002 May 27), 
which was accompanied by an M2 flar e but did not lead to a 
CME, was described by Ji et al.l J2003I) . The filament started 
to rise rapidly and developed a clear helical shape, as is often 
observed; however, the ascent was terminated at a projected 
height of « 80 Mm (Figs .[I]|2l. Such confined filament erup- 
tions are not uncommon (lRusi2003h . 
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Fig. 2. — Comparison of height, h(t), and velocity, u(t), of the flux rope 
apex in the simulation (solid lines; initial perturbation is dotted) with the 
corresponding values of the filament eruption in Fig. (data from Fig. 3 
of'ji et al. 12003]) overplotted as diamonds (the height data observed before 
18:04 UT were smoothed here, resulting in reduced velocity scatter). See 
text and Table fTI for the scaling of the dimensionless simulation variables 
(left axes) to the observed values (right axes). 



In order to model this event, we consider a kink-unstable 
configuration which is very similar to the case of an average 
flux rope twist of 4.9tt studied in detail in Paper I. The line 
current is reduced by about one third to enable a higher rise, 
but the average twist, $ = S.Ott here, is kept to reproduce the 
helical shape. (In the TD model, this fixes the minor radius 
to ^7 = 0.29, 8 percent larger than in the reference run in Pa- 
per I.) The sign of the line current is chosen to be positive 
to conform to the apparent positive (right-handed) helicity of 
the observed filament. Furthermore, we increase the numeri- 
cal diffusion and prevent the density from becoming negative, 
which permits us to follow the evolution of the system for a 
considerably longer time. Otherwise, the magnetic configura- 
tion and the numerical settings are the same as in Paper I. 

As in our previous simulations, the upwardly directed kink 
instability leads to the ascent and helical deformation of the 
flux rope as well as to the formation of current sheets (see 
Fig. 3 in Paper I). In Fig.^we compare the evolution of the he- 
lical shape of the flux rope with Transition Region and Coro- 
nal Explorer (TRACE) observations of the filament eruption. 
The evolution is remarkably similar Figure|2l shows that the 
principal features of the observed rise are also matched. After 
an exponential rise the flux rope comes to a stop at w 3.5 ho. 
A first deceleration occurs as the current density in the heli- 
cal current sheet above the apex begins to exceed the current 
density in the flux rope (f > 22). The subsequent upward push 
(f > 30) results from the reconnection outflow in the vertical 
current sheet below the rope. Finally, the rise is terminated by 
the onset of magnetic reconnection in the current sheet above 
the rope, which progressively cuts the rope field lines (f > 33). 
The reconnection outflows expand the top part of the rope in 
lateral direction, as seen both in observation and simulation. 

Using the scaling to dimensional values given in Tabled 
good quantitative agreement with the rise profile is obtained 
(Fig. 13, and the release of magnetic energy in this run of 
5 percent corresponds to 10^' erg, a reasonable value for a 
confined M2-class flare. 

The agreement between the observations of the event and 
our simulation confirms the long-held conjectures that the de- 
velopment of strongly helical axis shapes in the course of 
eruptions can be regarded as an indication of the kink instabil- 
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TABLE 1 
Parameter Settings 





Simulation Parameters 


Scaling Parameters 


Sect. 


<I>/7r 


b »y(2) 


L 




|Bo(/io)| W 










(Mm) (s) 


(G) (erg) 


3 


5.0 


0.29 0.83 


10 


23 11.5 


200 lO^i 


4 


-5.0 


0.33 1.54 


32 


115 111 


10-40 1031-32 



Note. — The expression ri(z) = -zt/lnBex(0, 0,z,0)/rfz is the 'decay in- 
dex' of the 'external' field (excluding the contribution by the ring current), L 
is the box size, and W is the released magnetic energy. The runs are equal 
in grid resolution in the central part of the box, A = 0.02, major rope ra- 
dius, R = 1 .83, and distance of the fictitious magnetic charges from the z axis, 
/ = 0.83. 
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Fig. 3. — Normalized initial magnetic field strength (thick lines) and Alfven 
velocity (thin lines) vs. height for the configurations described in Sect.l3l 
(original TD model; dashed lines) and Sect.|3 (modified TD model; solid 
lines). 



ity of a twisted flux rope and that the frequently observed heli- 
cal fine structures in erupting filaments and prominences out- 
line twisted fields. Furthermore, it shows that flux ropes with 
substantial twist can exist or be formed in the solar corona at 
the onset of, or prior to, eruptions. 

4. SIMULATION OF AN EJECTIVE ERUPTION (CME) 

The full eruption of the kink-unstable flux rope in the TD 
model is prevented by the strong overlying field, which is 
dominated by the line current. It is possible to obtain an erup- 
tive behaviour of the flux rope by removing the line current; 
however, such a modification leads to an infi nite number of 
field line turns at the surface of the flux rope jRoussev et alJ 
l2003h . In order to avoid this problem, we replaced the 
line current by a pair of subphotospheric dipoles (as used in 
|TCr6k & Kliem 2003). The position of the dipoles is chosen 
such that the field lines of the dipole pair passing through the 
flux rope match the curvature of the rope as closely as possi- 
ble. The resulting equilibrium yields a finite twist everywhere 
in the system, but the magnetic field overlying the flux rope 
now decreases significantly faster with height than in the orig- 
inal TD model (Fig.|3jl. By varying the dipole moments or the 
minor radius b, one can adjust the average twist within the 
flux rope. 

Choosing suitable dipole moments and b = 0.6, but oth- 
erwise the same parameters of the TD model as in Sect.|3j 
we first checked that the modified configuration relaxes to a 
nearby stable equilibrium for subcritical twist, <i> = 2.77r. Next 
a configuration with supercricital twist, $ = -S.Ott, is consid- 
ered, obtained by changing the minor radius to = 0.33 and 
reversing the dipole moments. The sign of the helicity corre- 
sponds to the 2001 May 15 event discussed below; it has no 
influence on the rise, h{t), of the flux rope apex. The numer- 
ical parameters of the simulation are the same as in Sect.|3l 
except for a considerably larger simulation box and a smaller 
level of numerical diffusion, which this system permitted. 

The helical kink instability also develops in the modified 
model. However, the flux rope now exhibits a much stronger 
expansion (Fig.E}, which is not slowed down. The initially 
exponential rise is followed by a rise with approximately con- 
stant and locally super-Alfvenic velocity, until the rope en- 
counters the top of the simulation box (at t w 80). The helical 
current sheet remains very weak on top of the flux rope apex 
and no significant amount of reconnection occurs here. On the 
other hand, the vertical current sheet now steepens in a large 
height range. Magnetic reconnection commences in this sheet 
at the beginning of the exponential phase and rises in tandem 
with the ascent of the flux rope, particularly closely during the 



exponential phase. Since the flux rope expands continuously 
during this phase (instead of being compressed by the upward 
reconnection outflow below it) and moves away from the fore- 
front of the outflow region afterwards, the ideal instability of 
the flux rope appears to be the driver of the closely coupled 
processes. Cusp-shaped field lines are formed throughout the 
evolution (Fig.|4} but most prominently in the late phase, in 
agreement with soft X-ray observations of eruptive flares. 

The full eruption of the flux rope in the modified TD model 
must be enabled by the weaker overlying field, since all other 
parameters are identical, or very close, to Sect.|3] 

In Fig.|5lwe compare the rise of the flux rope apex in the 
simulation with the rise of the apex of a well observed erup- 
tive prominence on 2001 May 15, which occurred in a spot- 
less region slightly behind the limb and led to a fast CME 
(peak velocity of leading edge « 1200 kms~' ) and a long- 
duration flare (X-ray class C4). As described bv'Marici c et alJ 
(2004), the eruptive prominence developed a helical shape, 
analogous to the middle panels in Fig. ^ with reversed hand- 
edness, and it showed the typical rise characteristics of a fast 
CME (initially exponential or e xponential-h ke rise, followed 
by approximately linear rise; Vrsnak "200?). For this com- 
parison, we first scaled the Alfven time such that the dura- 
tion of the exponential rise phase is matched, ta = 111 s, 
and shifted the time axis accordingly. Then we scaled the 
length unit such that the apex height at the point of peak 
acceleration in the simulation it = 21) equals the height of 
the prominence at the resulting observation time (see bottom 
axis), i.e., /zq = 115 Mm. This fixes the scaling of the velocity 
and acceleration amplitudes. Apart from a somewhat more 
gradual decrease of the observed acceleration after the peak, 
excellent qualitative and quantitative agreement is obtained, 
demonstrating (as in Sect. |3} that the kink instability yields 
the growth rate required by the observed rise profile for the 
twist indicated by the observed helical shape. The slight dif- 
ference in the late acceleration profile may have many origins, 
for example, a different height profile of the field strength, or 
a simultaneous expansion of the overlying field enforced by 
photospheric flows in the observed event. 

The simulation shows a strong magnetic energy release, 
25 percent of the initial value, which agrees with the mag- 
nitude observed in ejective solar eruptions jForbesI [ZOOCt 
Emslie et al. 2004). The considered event may have released 
magnetic energy of order ^ (lO-^'-lO-'^) erg (the low X-ray 
class resulted from footpoint occultation). Using the above 
scaling for /zq, this energy release is reproduced for Bq ^ 
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Fig. 4. — Magnetic field fines of the kink-unstable modified TD model at ? = (left), t = 30 (center), and ; = 43 (right). The magnetogram, B,(x,y,0,t), is 
included. Field lines started at a circle of radius h/3 in the bottom plane show the core of the flux rope. Additional green field fines, also with identical start 
points in all panels, indicate the formation of "post flare loops" with a cusp by reconnection. The hyperbolic point of the field at the z-axis (magnetic X-point) 
lies a.tz~ 0.2, 0.6, and 1.1, respectively. 
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Fig. 5. — Comparison of the simulation in Sect.|4]with the CME on 2001 
May 15 in the same format as in Fig.|2] including the acceleration a(t). Dia- 
monds show the ape x motion of the erupti ng prominence/the CME core (data 
are from Fig. 6a, c of lMaricic et alJ I2004r '). 



of solar eruptions, identify a mech anism for the p rocess, con- 
firming the original suggestion bv^ Sakurail (Il976h . They also 
demonstrate the importance of the height dependence of the 
overlying field for the evolution of the instability into a CME, 
while Amari et al. (2003b) found the amount of magnetic he- 
licity to be essential. Since lAmari et al. built up the helicity 
by rotating the main photospheric polarities, which simulta- 
neously expands the overlying field (Torok & Kliem 2003J, 
both results are fully consistent with each other 



5. CONCLUSIONS 

Our MHD simulations of the kink instability of a coronal 
magnetic flux rope reproduce essential properties — an ini- 
tially exponential rise with the rapid development of a heli- 
cal shape — of two well observed solar eruptions, one of them 
confined, the other ejective. The subsequent approximately 
linear rise of the ejective eruption is reproduced as well. Since 
these f eatures are ch aracteristic properties of many solar erup- 
tions ( Vrsnakl2001l) . we regard the kink instability of coronal 
magnetic flux ropes as the initiation mechanism and initial 
driver of many such events. A sufficiently steep decrease of 
the magnetic field with height above the flux rope permits the 
process to evolve into a CME. 



(10^0) G, consistent with expected averages of the coronal 
field strength over the large length scales involved. 

The simulation could also be scaled to a filament eruption 
that was associated with an X-class flare and a ve ry fast CME 
(on 2004 November 10; see 'Williams et al.l2005h . 

A line-tie d flux rope w as found to erupt in a few previous 
simulations ("Amari et al. 2000, 2003a b). The present simu- 
lations, through their agreement with characteristic properties 
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